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ABSTRACT 

We investigate the luminosity functions of Lyman-break galaxies (LBG) at z ~ 4 and 5 
based on the optical imaging data obtained in the Subaru Deep Field (SDF) Project, a program 
conducted by Subaru Observatory to carry out a deep and wide survey of distant galaxies. Three 
samples of LBGs in a contiguous 875 arcmin^ area are constructed. One consists of LBGs at 
z ~ 4 down to i' — 26.85 selected with the B — Rvs R—i' diagram (BRi'-LBGs). The other two 
consist of LBGs at z ^ 5 down to z' = 26.05 selected with two kinds of two-color diagrams: V — i' 
vs i'-z' (Vi'z'-LBGs) and R-i' vs i'-z' (ffi'z'-LBGs). The number detected is 3, 808 for BRi'- 
LBGs, 539 for yi'z'-LBGs, and 240 for i?i'z'-LBGs. The adopted selection criteria are proved to 
be fairly reliable by the spectroscopic observation of 63 LBG candidates, among which only 2 are 
found to be foreground objects. We estimate the fraction of contamination and the completeness 
for these three samples by Monte Carlo simulations, and derive the luminosity functions of the 
LBGs at rest-frame ultraviolet wavelengths down to Afyv — —19.2 at z '--^ 4 and A/uv = —20.3 
at z 5. We find clear evolution of the luminosity function over the redshift range of < z < 6, 
which is accounted for by a sole change in the characteristic magnitude, M*. The cosmic star 
formation rate (SFR) density at z ~ 4 and z 5 is measured from the luminosity functions. 
The measurements of the integrated SFR density at these redshifts are largely improved, since 
the luminosity functions are derived down to very faint magnitudes. We examine the evolution 
of the cosmic SFR density and its luminosity dependence over < z < 6. The SFR density 
contributed from brighter galaxies is found to change more drastically with cosmic time. The 
contribution from galaxies brighter than M*^^ — 0.5 has a sharp peak around z = 3 - 4, while 
that from galaxies fainter than M*^^ — 0.5 evolves relatively mildly with a broad peak at earlier 
epoch. Combining the observed SFR density with the standard Cold Dark Matter model, we 
compute the cosmic SFR per unit baryon mass in dark haloes, i.e., the specific SFR. The specific 
SFR is found to scale with redshift as (1 -I- z)^ up to z ^ 4, implying that the efficiency of star 
formation is on average higher at higher redshift in proportion to the cooling rate within dark 
haloes, while this is not simply the case at z > 4. 

Subject headings: cosmology: observations — galaxies: evolution — galaxies: high-redshift — galaxies: 
luminosity function, mass function 



1 



1. INTRODUCTION 

When and how galaxies formed is one of the 
primary questions in astronomy today. Obser- 
vations of young galaxies at high redshifts are a 
straightforward approach to this problem. Over 
the past decade, it has become feasible to under- 
take large surveys of galaxies at high redshifts. 
What made this possible is the progress in observ- 
ing technology including telescopes and detectors, 
and also sophistication of selection methods to lo- 
cate high-redshift galaxies. Probably the most ef- 
ficient method is the Lyman-break technique. It is 
a simple photometric technique based on the con- 
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tinuum features in rest-frame ultraviolet spectra 
redshifted into optical bandpasses, and requires 
only optical imaging in a few bandpasses. Pio- 
neered by Guhathakurta et al. (1990), this method 
has been successfully used to find many young, 
star-forming galaxies beyond ^; ~ 2 (e.g., Steidel 
& Hamilton 1992; Steidel ct al. 1999; Lehnert & 
Bremer 2003; Iwata et al. 2003; Ouchi et al. 2004; 
Dickinson et al. 2004; Sawicki & Thompson 2006). 
Steidel et al. (2003) made spectroscopic observa- 
tion for about 1000 photometrically selected z ^ 3 
galaxies and verified the usefulness of this method. 
High-redshift galaxies selected by this method are 
called Lyman-break galaxies (LBGs). 

Since their first discovery in the 1990s, various 
properties of LBGs have been extensively studied. 
One of the most fundamental measurements is the 
luminosity function (LF) at rest-frame ultraviolet 
wavelengths (i.e., the number density of galaxies 
as a function of ultraviolet luminosity). Besides 
providing information on the number density of 
galaxies, the LF can be used to probe the star for- 
mation activity in the universe, because ultravio- 
let luminosity is sensitive to star formation. Thus 
the ultraviolet luminosity density derived by in- 
tegrating the LF is related to the star formation 
rate (SFR) density in the universe. By obtaining 
the LF at various redshifts and examining its evo- 
lution, one can gain insights into the formation 
history of galaxies and the star formation history 
in the universe. 

Steidel ct al. (1999) derived the LF at z - 3 
from their large LBG sample, finding that their 
data are well fitted by a Schechter function with a 
low-luminosity slope a of —1.6 down to i ~ O.IL* . 
They extended their LBG search to z ~ 4, to de- 
tect no significant evolution of the LF at bright 
magnitudes (M < —21.1) from z ^ 4 to z ^ 3. 
More recent observations have also derived the LF 
at z ^ 4 and z ~ 5 (Iwata ct al. 2003; Ouchi et 
al. 2004; Gabasch ct al. 2004; Sawicki & Thomp- 
son 2006). The frontier of LBG searches now ex- 
tends beyond z ~ 5. (Yan et al. 2003; Stanway et 
al. 2003; Bunker et al. 2004; Bouwens et al. 2004, 
2005; Shimasaku et al. 2005). 

However, it is not easy to construct a large sam- 
ple of LBGs beyond z 4 covering a wide range 
of absolute magnitude because of their apparent 
faintness and low surface number density. In the 
first place, most of the studies to date are limited 
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to bright LBGs. Consequently, the shape of the 
faint end of the LF has not been well determined. 
Thus, the estimation of the cosmic SFR density 
has a considerably large uncertainty (about a fac- 
tor of 3 - 10) due to the long extrapolation of the 
observed LFs to faint magnitudes (see Ouchi et al. 
2004). Exploring LBGs down to faint luminosities 
is necessary in order to determine the overall shape 
of the LF, and to measure the cosmic SFR density 
accurately. Although some surveys arc extremely 
deep, they are restricted to very small areas (e.g. 

4 arcmin^ for the Hubble Deep Field, and ~ 10 
arcmin^ for the Hubble Ultra Deep Field). Sur- 
veys based on such a small area probably suffer 
from cosmic variance, i.e., inhomogeneities in the 
spatial distribution of LBGs. Although deep and 
wide surveys of LBGs at z ~ 4 have been made by 
Gabasch et al. (2004) and Sawicki & Thompson 
(2006) very recently, large samples for z > 4 arc 
still very limited. It is therefore crucial to con- 
struct a new LBG sample from a survey of a sim- 
ilar depth and width and derive the LF and the 
cosmic SFR independently. 

In this paper, we present a detailed study of 
LBGs at z ~ 4 - 5 based on the deep and wide-field 
images obtained in the Subaru Deep Field (SDF) 
Project (Kashikawa et al. 2004). The SDF Project 
is a program conducted by Subaru Observatory to 
carry out a deep galaxy survey over a blank field 
as large as ~ 0.2 square degree. Exploiting these 
unique data, we make the largest samples of LBGs 
at 2: ~ 4 and 5, and derive the LFs down to very 
faint magnitudes; Muv = —19.2 at 2; ~ 4 and 
Muv = -20.3 at 2 ~ 5. This extends the LBG 
search by Ouchi et al. (2004) based on shallower 
and smaller area data of the same field some ~ 
0.5 mag further and 1.5 times wider, leading to a 
~ 3 times increase in number of LBGs detected. 
These samples enable us to examine the behavior 
of the LF more accurately over a wide magnitude 
range and obtain more reliable measurements of 
the cosmic SFR density in the early universe. The 
LFs of LBGs at z ^ 4 and 5 we derive arc found 
to differ from those at the same redshift ranges 
obtained in some of the previous surveys. 

Recently, the LF at ultraviolet wavelengths of 
present-day galaxies {z ~ 0) has been derived very 
accurately from a large galaxy survey by Galaxy 
Evolution Explorer (GALEX) (Wyder et al. 2005). 
Arnouts et al. (2005) have derived the LFs at 



0.2 < z < 1.2 also based on GALEX observa- 
tions. Measurements at z ~ 6 are now avail- 
able (Bouwens et al. 2005; Shimasaku et al. 2005). 
Combining the results at z ^ 4 and 5 from this 
work with those at lower redshifts and z ~ 6 from 
the literature, we investigate the evolutionary his- 
tory of galaxies in terms of the star formation ac- 
tivity over < z < 6. 

The outline of this paper is as follows. Sec- 
tion 2 describes the data of the SDF Project used 
in this study. Section 3 describes the selection of 
LBGs at 2; ~ 4 and 5 from the photometric data. 
The contamination by interlopers and the com- 
pleteness of the samples are also estimated. In 
Section 4, the rest-frame UV luminosity functions 
of LBGs at z ^ 4 and 5 are derived, and com- 
pared with those by other authors. In Section 5 
we discuss the evolution of the cosmic star forma- 
tion activity over the redshift range < z < 6. 
The evolution of the luminosity function and the 
cosmic star formation rate density is examined. 
Section 6 summarizes this study. 

Throughout this paper, the photometric system 
is based on AB magnitude (Oke & Gunn 1983). 
The cosmology adopted is a flat universe with 
= 0.3, r^A = 0.7, and a Hubble constant of 
Hq = 70 km Mpc~^. These values are con- 
sistent with those obtained from the latest CMB 
observations (Spergel et al. 2003). 

2. DATA 

The data used in this study are obtained in 
the SDF Project (Kashikawa et al. 2004). The 
program consists of very deep multi-band opti- 
cal imaging, NIR imaging for smaller portions of 
the field, and follow-up optical spectroscopy. This 
study is based on the optical imaging data to- 
gether with additional information from the spec- 
troscopic data. Full details of the optical imaging 
observations, data reduction, and object detection 
and photometry are described in Kashikawa et al. 
(2004). 

2.1. Imaging Data 

The optical imaging observations in the SDF 
Project were carried out with the prime-focus 
camera (Suprime-Cam: Miyazaki et al. 2002) 
mounted on the 8.2 m Subaru Telescope (lye et 
al. 2004) in 2002 and 2003. The imaging was 
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made for a single field of view of the Suprime-Cam 
(34' X 27') toward the Subaru Deep Field (SDF: 
Maihara et al. 2001), centered on [13''24™38/9, 
+27°29'25."9 (J20GG)]. The scale of the Suprime- 
Cam is 0."202 pixel" ^. The images were taken 
in five standard broad-bands filters: B, V, R, i', 
and z' , and two narrow-bands filters: NB816 and 
NB921. The central wavelengths of the broad- 
bands are 4458 A, 5478 A, 6533 A, 7684 A, and 
9037 A for B, V, R, i', and z', respectively. The 
SDF was previously imaged during the commis- 
sioning runs of the Suprime-Cam in 2001 in five 
bandpasses, B, V , R, i', and z', with 1 3 hours 
exposure times. The work by Ouchi et al. (2004) 
is based on these data. They were co-added in 
constructing the final images. 

All individual CCD data with good quality 
were reduced and combined to make a single im- 
age for each band using a pipeline software pack- 
age (Ouchi 2003) whose core programs were taken 
from IRAF and mosaic-CCD data reduction soft- 
ware (Yagi et al. 2002). The combined images for 
all seven bands were aligned and smoothed with 
Gaussian kernels so that all have the same see- 
ing size, a PSF FWHM of 0."98. The efficiency 
and reliability of object detection and photome- 
try in regions near the edges of the images are 
significantly lower on account of low S/N ratios 
caused by dithering. Regions around very bright 
stars are also degraded due to bright haloes and 
saturation trails. We carefully defined these low- 
quality regions, and did not use objects located 
in these regions. The effective area of the final 
images after removal of the low-quality regions is 
875 arcmin^. Photometric calibration for the im- 
ages was made using standard stars taken during 
the observations. The total exposure time reaches 
~ 10 hours for each band, and the limiting magni- 
tude (3fT within a 2" diameter aperture) is 28.45 
(B), 27.74 {V), 27.80 (R), 27.43 (i'), 26.62 {z'), 
26.63 (NB816), and 26.54 (NB921). 

Object detection and photometry were per- 
formed using SExtractor version 2.1.6 (Bertin & 
Arnouts 1996). A collection of at least 5 contigu- 
ous pixels above the 2a sky noise was identified as 
an object. The object detection was made for all 
seven images independently. For each of the ob- 
jects detected in a given band image, photometry 
was made in all the images at exactly the same 
position as in the detection-band image. Thus, 



seven catalogs were constructed for respective de- 
tection bands. For the present study, we limit the 
catalogs to objects whose detection-band magni- 
tudes are brighter than the 5(7 limiting magnitude 
(5(7 sky noise within a 2" diameter aperture) of 
that band, in order to provide a reasonable level 
of photometric completeness. We use magnitudes 
within a 2" diameter aperture to derive the colors 
of objects, and adopt MAG_AUTO in SExtractor 
for total magnitudes. The magnitudes of objects 
were corrected for a small amount of foreground 
Galactic extinction toward the SDF using the dust 
map of Schlegel et al. (1998). Since the variation 
of E{B—V) over the SDF of 875 arcmin^ is at most 
0.007, we adopt a single value, E{B -V) = 0.017, 
for correction, which is the value at the center of 
the SDF. The value of E{B -V) = 0.017 corre- 
sponds to an extinction of Ab = 0.07, Ay = 0.05, 
Ar = 0.04, A^^ = 0.03, A,, = 0.02, Anb816 = 0.03, 
and ^NB92i = 0.02. 

2.2. Spectroscopic Data 

We carried out spectroscopic follow-up obser- 
vations for some objects in our catalogs with the 
Subaru Faint Object Camera and Spectrograph 
(FOCAS: Kashikawa et al. 2002) on the Subaru 
Telescope during 2002 - 2004 and DEIMOS (Faber 
et al. 2003) on the Keck II Telescope in 2004. The 
spectroscopic observations were made in the multi- 
object slit mode. For the FOCAS observations, 
we used a 300 line mm~^ grating with an 058 
order-cut filter and an Echelle with a z' filter. Slit 
widths on the masks were either 0."8 or 0."6 for 
the 300 line configuration, resulting in a spectral 
resolution of 9.5 A and 7.1 A at 8150 A, respec- 
tively, and 0."8 for the Echelle. The integration 
times of respective masks were 7000 - 9000 sec- 
onds. Flux calibration was made with spectra of 
the spectroscopic standard stars Hz 44 and Feige 
34. The data were reduced in a standard man- 
ner. For the observations with DEIMOS, we used 
a 830 line mm~^ grating with a GG495 order-cut 
filter. Slit widths were l."0, which gave a spec- 
tral resolution of 3.97 A. The integration times 
were 6300 - 19800 seconds. Spectra of the stan- 
dard stars BD-|-28d4211 and Feige 110 were taken 
for flux calibration. The data were reduced with 
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Fig. 1. — The magnitude distribution of the spec- 
troscopic objects with .Zgpec < 5.5 in our catalogs. 



Fig. 2. — The redshift distribution of the spec- 
troscopic objects with .Zgpec < 5.5 in our catalogs. 



the spec2d pipeline ^ for the reduction of DEEP2 
DEIMOS data. In addition, we also use spectra 
of objects in the SDF which were taken during 

the guaranteed time observations of the FOCAS 
in 2001 (Kashikawa ct al. 2003; Ouchi et al. 2004). 

In Figures 1 and 2, we present the distribution 
of the z' magnitude of the spectroscopic objects 
and their redshift distribution, respectively. Only 
objects with Zgpec < 5.5 are included here. 

3. LYMAN-BREAK GALAXY SAMPLES 

AT z ~ 4 AND 5 

3.1. Selection of Lyman-break Galaxies 

We demonstrate the principle of the Lyman- 
break color selection technique to locate galaxies 
at z ^ 4 and z ^ 5 in Figure 3. For galaxies at 
z ~ 4, the break in the spectrum at rest-frame ul- 
traviolet wavelength enters the B band and the flat 
continuum longward of Lyman-a shifts to wave- 
lengths longer than the R band, and these galax- 
ies are identified by their red B — R color and blue 
R—i' color. We select LBGs at z 4 by these two 
colors (hereafter called BRi'-LBGs). For galaxies 
at z ~ 5, the spectral break enters into the V band 
and the flat continuum is sampled at wavelengths 
longer than the i' band. At z ~ 5, the attenuation 



^The spec2d pipeline was developed at UC Berkeley with 
support from NSF grant AST-0071048. 



due to the Lyman-a forest is so strong that the 
flux in the R band is severely suppressed, as well 
as the flux in the V band. Thus, galaxies at z ~ 5 
arc identified by their red V — i' or R—i' color and 
blue i' — z' color. We make two samples for LBGs 
at z ~ 5: one selected by V—i' and i' — z' colors, 
and the other selected by R~i' and i' — z' colors 
(hereafter called Vi' z'-LQGs and Ri'z'-LQGs, re- 
spectively). 

Figures 4-6 illustrate the predicted posi- 
tions of high-rcdshift galaxies and foreground ob- 
jects (lowcr-redshift galaxies and Galactic stars) 
in three two-color diagrams: B — R vs R — i' , V — i' 
vs i' — z', and R — i' vs i' — z' . The solid lines in 
these figures indicate the tracks for model spec- 
tra of young star-forming galaxies for the redshift 
range z > 3. The model spectrum was constructed 
using the stellar population synthesis code devel- 
oped by Kodama & Arimoto (1997). For model 
parameters, an age of 0.1 Gyr, a Salpeter ini- 
tial mass function, and a star-formation timescale 
of 5 Gyr were adopted, and reddening of E{B — 
V) = 0.16 was applied using the dust extinction 
formula for starburst galaxies by Calzetti et al. 
(2000). These values reproduce the average rest- 
frame ultraviolet-optical spectral energy distribu- 
tion of LBGs observed at z ~ 3 (Papovich et al. 
2001). We also show the tracks for model spec- 
tra with reddening of E{B — V) = and 0.3 for 
reference. The absorption due to the intergalac- 
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Fig. 3. — Upper panel: The B, R, and i' baud- 
passes overplotted on the spectrum of a generic 
2: = 4 galaxy (thick hne), illustrating the utility of 
color selection technique using these three band- 
passes for locating z ~ 4 galaxies. Lower panel: 
The V, R, i', and z' bandpasses overplotted on the 
spectrum of a generic z = 5 galaxy (thick line). 
The band sets of Vi'z' and Ri'z' work well to iso- 
late 2; ~ 5 galaxies. 

tic medium was applied following the prescription 
by Madau (1995). The dotted, dashed, and dot- 
dashed lines delineate the tracks for model spec- 
tra of local elliptical, spiral, and irregular galax- 
ies, respectively, redshifted from z = to z = 3. 
These spectra were also constructed using the stel- 
lar population synthesis code by Kodama & Ari- 
moto (1997), and redshifted without evolution. 
The asterisks denote the 175 Galactic stars of 
various types whose spectra are given by Gunn 
& Stryker (1983). The colors are calculated by 
convolving constructed model spectra of galaxies 
or given spectra of stars with the response func- 
tions of the Suprime-Cam filters. It can be seen 
that the model young star-forming galaxy moves 
nearly vertically from the origin in these two-color 
diagrams at z > 3, namely it becomes redder 
in B — R, V — i', and R—i' colors, while keeping 
R—i', i' — z', and i' — z'colois blue, as redshift in- 
creases, and that it shifts into a portion of the two- 
color space which is unpopulated by foreground 



objects. At z > 5, the R~i' color does not be- 
come redder very much while i' — z' color becomes 
redder in the R—i' vs i' — z' diagram. These fig- 
ures imply that LBGs at z ^ 4 and those at z ~ 5 
can be well isolated from foreground objects in 
the B — R vs R—i' diagram, and in the V—i' vs 
i' — z' and R—i' vs i' — z' diagrams, respectively. 
The most critical contamination is caused by el- 
liptical galaxies at z ~ 0.5, which come close to 
high-redshift galaxies in these diagrams due to the 
4000 A break on their spectra. 
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Fig. 4. — B — R vs R—i' diagram where the pre- 
dicted colors of model galaxies and stars are illus- 
trated. The three solid lines indicate the tracks 
for model spectra of young star-forming galaxies 
with reddening of E{B - y) = 0, 0.16, and 0.3 
(from left to right). 

The redshift range is from z = 3 to higher red- 
shifts, and the circles on the track mark the red- 
shift interval of 0.1. The dotted, dashed, and dot- 
dashed lines delineate the tracks for model spectra 
of local elliptical, spiral, and irregular galaxies, re- 
spectively, redshifted from z = to 3 without evo- 
lution. The circles on each track mark z = 0, 1, 
and 2. The asterisks represent the colors of 175 
Galactic stars given by Gunn & Stryker (1983). 
The thick line indicates the boundary which we 
adopt for the selection of BRi'-LQGs. 
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Fig. 5. — V — i' vs i' — z' diagram where the pre- 
dicted colors of model galaxies and stars are il- 
lustrated. Lines and symbols arc the same as in 
Figure 4. The thick line indicates the boundary 
which we adopt for the selection of Vi'z'-LBGs. 



Fig. 6. — R—i' vs i' — z' diagram where the pre- 
dicted colors of model galaxies and stars are il- 
lustrated. Lines and symbols arc the same as in 
Figure 4. The thick line indicates the boundary 
which we adopt for the selection of Ri'z'-LBGs. 



In Figures 7-9, we show the distribution of 
all the objects contained in the catalogs in the 
two-color diagrams. When the magnitude of an 
object in a non-detection band is fainter than the 
1(7 magnitude of the band, the la magnitude is as- 
signed to the object. Objects with spectroscopic 
redshifts are shown with colored symbols in these 
figures, where diS^erent colors mean diflferent red- 
shift bins. It is found that the colors of spectro- 
scopically identified objects match fairly well with 
those of model galaxies with corresponding red- 
shifts in Figures 4-6. 

We adopt the same color criteria for BRi' -LBGs 
and Vi'z'-LBGs as those of Ouchi et al. (2004). 
For Ri' z'-LBGs, we visually fine-tune their color 
criteria based on the increased redshift informa- 
tion. Specifically, we set the selection criteria for 
BRi'-LBGs as: 



B-R> 1.2, 
R-i' < 0.7, 

B -R>1.6{R-i') + 1.9, 



(1) 



for yi'^'-LBGs as: 

V-i' > 1.2, 
i' -z' < 0.7, 
V-i' > 1.8{i' 

B > 3a mag, 
and for Ri'z'-LBGs as: 



z') + 1.7, 



R-i' > 1.0, 
i' ~z' < 0.7, 

R-i' > 1.2{i' - z') + 0.9, 
B, V >3a mag. 



(2a) 



(2b) 



(3a) 



(3b) 



The boundaries on the two-color diagrams defined 
by these color criteria are outlined with the thick 
orange lines in Figures 4 6 and 7 -9. For the 
bands which are placed entirely shortward of the 
Lyman break, we additionally require that LBGs 
should be non-detected ((2b) and (3b)), which is 
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Fig. 7. — B — Rvs R—i' diagram for all the objects 
contained in the i'-detection catalogs (106,025 ob- 
jects down to i' < 26.85). When the B and/or 
R magnitude of an object is fainter than the la 
magnitude of the band, the la magnitude is as- 
signed to the object. The horizontal sequence 
along B — R = 0.62 shows objects which are fainter 
than the la magnitude in both B and R. 
The colored symbols show objects with spectro- 
scopic redshifts, where cyan, green, violet, and red 
represent objects in the range z < 3.0, 3.0 < z < 
3.7, 3.7 < 2; < 4.0, and 4.0 < z < 4.5, respec- 
tively. The thick orange line indicates the bound- 
ary which we adopt for the selection of BRi'- 
LBGs. 



Fig. 8. — V — i' vs i' — z' diagram for all the objects 
contained in the z'-detection catalogs (76,351 ob- 
jects down to z' < 26.05). When the V and/or 
i' magnitude of an object is fainter than the la 
magnitude of the band, the la magnitude is as- 
signed to the object. The horizontal sequence 
along V — i' ~ 0.29 shows objects which are fainter 
than the la magnitude in both V and i' . 
The colored symbols show objects with spectro- 
scopic redshifts, where cyan, green, violet, and red 
represent objects in the range z < 3.0, 3.0 < z < 
4.4, 4.4 < z < 4.7, and 4.7 < z < 5.5, respec- 
tively. The thick orange line indicates the bound- 
ary which we adopt for the selection of Vi'z'- 
LBGs. 



expected to work effectively to reduce contamina- 
tion from foreground galaxies. 

We apply these selection criteria to the cata- 
logs. We use the i'-detection catalogs for BRi'- 
LBGs, which contain 106, 025 objects to a limit 
of i' — 26.85 (total magnitude corrected for the 
Galactic extinction), and the z'-detection cata- 
logs for T^«'z'-LBGs and Ri' z'-LBGs, which con- 
tain 76,351 objects to z' — 26.05 (total magni- 
tude corrected for the Galactic extinction). The 
number of objects which satisfy the selection cri- 
teria, i.e., the number of LBG candidates, is 3, 808 
for BRi'-LBGs, 539 for Vi'z'-LBGs, and 240 for 



Ri' z'-LBGs ^. Tabic 1 summarizes the samples of 
the LBG candidates. 

The selection criteria are shown by our spec- 
troscopy to be fairly reliable. We spectroscopi- 
cally observed 105 objects in the BRi'-LBG sam- 
ple, and identified 42 objects. Among the 42, only 
one is an interloper at z = 3.29, with the remain- 
ing 41 being at z > 3.7. Similarly, 32 objects from 

^For this study, we limit the catalogs to the 5ct Umiting mag- 
nitude. Kashikawa et al. (2006) also use our LBG samples 
for their study on clustering properties of LBGs but down 
to the So- limiting magnitude. The number of LBG candi- 
dates with magnitudes brighter than the 3cr limiting mag- 
nitude is 4, 543 for Bffi'-LBGs, 831 for Vi'z'-LBGs, and 
386 for ffi'z'-LBGs. This difference causes no significant 
effect on luminosity functions except for the faintest bins. 
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Fig. 9. — R—i' vs i' — z' diagram for all the objects 
contained in the z'-detection catalogs (76,351 ob- 
jects down to z' < 26.05). When the R and/or 
i' magnitude of an object is fainter than the la 
magnitude of the band, the la magnitude is as- 
signed to the object. The horizontal sequence 
along R—i' — 0.35 shows objects which are fainter 
than the la magnitude in both R and i' . 
The colored symbols show objects with spectro- 
scopic redshifts, where cyan, green, violet, and red 
represent objects in the range z < 3.0, 3.0 < z < 
4.7, 4.7 < z < 4.9, and 4.9 < z < 5.5, respec- 
tively. The thick orange line indicates the bound- 
ary which we adopt for the selection of Ri'z'- 
LBGs. 



lowing simulation, however, we infer that the frac- 
tion of such interlopers is not so high in our LBG 
sample. First, we add noise to the colors of all ob- 
jects outside of the BRi' selection criteria in our 
photometric catalog according to their measured 
magnitudes. We then count the objects which 
now meet the selection criteria due to the mod- 
ified colors, to find that they are less than 10% of 
the original LBG candidates even for the faintest 
magnitude bin where photometric errors are the 
largest. Note that this simulation gives an upper 
limit to the fraction of interlopers, since the mea- 
sured colors already include photometric errors. 
In Figure 10, we present the apparent magnitude 
distribution of the spectroscopically identified ob- 
jects in each LBG sample. Although the targets 
of spectroscopy is biased toward bright objects, 
we emphasize here that the spectroscopic samples 
also include faint objects and there are no interlop- 
ers even among them. We cannot, however, rule 
out the possibility that there may be interlopers 
for which no redshift have been secured. Figure 11 
shows the redshift distribution of the spectroscopic 
samples. The average redshift is < z >spcc= 4.1 
for BRi'-LBGa, < z >spec= 4.8 for Vi'z'-LBGs, 
and < z >spcc= 4.9 for i?i'z'-LBGs. Examples of 
spectra of z ~ 4 LBGs and z ~ 5 LBGs are shown 
in Figures 12 and 13, respectively. 



the V^j'z'-LBG sample were spectroscopically ob- 
served, and 23 were identified. All but one of the 
23 objects are at z > 4.4; the remaining one ap- 
pears to be a Galactic star. In the Ri'z'-LBG sam- 
ple, 12 objects were spectroscopically observed, 
and 9 were identified. All are found to be at 
z > 4.5. Possible reasons for the low detection 
rate for BRi'-LBGs are that the Lya emission of 
z ~ 4 LBGs is on average weaker than that of 
z 5 LBGs and that the Lya line at z ^ 4 falls in 
the wavelength range where the sensitivity of our 
spectroscopy is relatively low. The failed sample 
may also include interlopers between z ^ 2 and 3 
whose strong line features do not enter the wave- 
length coverage of our spectrograph. From the fol- 
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Table 1 
Samples of LBG candidates 



SMiuplc Najiu" 


Nimiljcr of CVuulidatoh 


DcttX'lioii Band 


Mag- 


iitu(l( 


riaiigc* 


BRi'-LBG 


3,808 


i' 


i' = 


22.85 


- 26.85 


Vi'z'-LBG 


539 


z' 


z' = 


23.55 


-26.05 


Ri'z'-LBG 


240 


z' 


z' = 


24.05 


- 26.05 



* Total magnitudes corrected for the Galactic extinction. 



15 



I 10 



13 





15 



% 10 



-Vi'z'-LBGs 




15 



0) 



3 




-ffi'z'-LBGs 



23 



24 25 
m 



_L 



26 27 



.£3 



13 



10 

8 

6 

4 

2 


10 



^5 



'-LBGs 



Wz'-LBGs 




i?i 'z'-LBGs 



3 3.5 -1 4.5 5 5.5 6 
Kedshift 



Fig. 10. — The magnitude distribution of the spec- 
troscopically identified objects in each LBG sam- 
ple. Dotted histograms show the magnitude dis- 
tribution of aU the objects in each LBG sample 
normalized to the number of the spectroscopically 
identified objects. 



Fig. 11. — The redshift distribution of the spectro- 
scopically identified objects in each LBG sample. 
Dotted lines show the redshift distribution func- 
tions from the Monte Carlo simulation (see §3.3) 
normalized to the number of the spectroscopically 
identified objects. 
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Fig. 12. — Examples of spectra of Lyman-break 
galaxies at 2 ~ 4 with good quality that have 
prominent spectral features such as a strong Lya 
emission line and interstellar lines. 
The lowest panel shows a relative night-sky spec- 
trum. 

3.2. Number Counts 

The observed number counts of the LBG can- 
didates as a function of apparent magnitude are 
shown in Figures 14 - 16, along with those mea- 
sured by other authors who selected similar LBGs. 
None of the data is corrected for contamination or 
incompleteness. Previous LBG surveys included 
in the figures are summarized in Table 2. Ouchi 
et al. (2004) used the same photometric systems to 
make samples of LBGs at ^; ~ 4 over a 543 arcmin^ 



Fig. 13. — Same as Figure 12, but for z ^ 5. 

area and LBGs at z ^ 5 over a 616 arcmin^ area of 
the SDF. Stcidcl ct al. (1999) carried out a search 
for LBGs at z 4 using G — 91, 91—/ colors in 
a total of 828 arcmin^ area of 10 separate fields. 
Sawicki & Thompson (2006) undertake a survey of 
LBGs at z ~ 4 with the same photometric systems 
and selection criteria that Steidel et al. (1999) use 
but to a deeper magniude in the Keck Deep Fields 
which cover a total area of 169 arcmin^ and consist 
of five fields. Iwata et al. (2003) found LBGs at 
z ~ 5 using V — Ic, Ic — z' colors in a 575 arcmin^ 
area including the Hubble Deep Field North. Ca- 
pak et al. (2004) obtained number counts of LBGs 
at z ~ 4 selected hy B — R, R—Ic colors and LBGs 
at z ~ 5 selected hyV — Ic,Ic — z' colors from a 720 
arcmin^ area centered on the Hubble Deep Field 
North. Our limiting magnitudes are among the 
faintest ones in the LBG samples for z ^ 4 and 5 
constructed to date. 

In general, number counts measurements of 
LBGs are dependent on detection completeness, 
photometric systems used, and selection crite- 
ria. It is found that our counts match very well 
with those of Ouchi et al. (2004) for all the sam- 
ples except for brightest bins. While we use 
MAG_AUTO in SExtractor as total magnitudes of 
objects, Ouchi et al. (2004) calculate total magni- 
tudes from aperture magnitudes through a fixed 
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Table 2 

lbg surveys in the literature 



Sample Area [arcmin^] JTiUmit ^ Selection Fig. in this paper Ref. 



z 


~ 4 


875 


i' 




26.85 


3,808 


B-R, R-i' 


Fie; 


14, 19, 21 


This study 


z 


~ 5 


875 


z' 




26.05 


539 


V-i', i'-z' 


Fig. 


15, 20, 21 


This study 


z 


~ 5 


875 


z' 





26.05 


240 


R-i', i'-z' 


Fig. 


16, 20, 21 


This study 


z 


~ 3 


1046 


m 





25.5 


1,270 


Un-G, G-fH 


Fig. 


21 


Steidel et al. (1999) 


z 


^4 


828 


I 




25.0 


207 


G-m, fR-/ 


Fig. 


14, 19 


Stcidcl ct al. (1999) 


z 


~4 


543 


i' 




26.3 


1,438 


B-R, R-i' 


Fig. 


14, 19 


Ouchi et al. (2004) 


z 


^4 


720 


Ic 




25.5 


? 


B — R, R — /c 


Fig. 


14 


Capak et al. (2004) 


z 


~4 


169 


91 




27.0 


427 


G-m, m-i 


Fig. 


19 


Sawicki & Thompson (2006) 


z 


~ 5 


575 


Ic 




25.95 


305 


v-ic, ic-z' 


Fig. 


15, 20 


Iwata et al. (2003) 


z 


~ 5 


616 


z' 




25.8 


246 


v-i', i'-z' 


Fig. 


15, 20 


Ouchi et al. (2004) 


z 


~ 5 


720 


Ic 




25.5 


? 


V-h, h~z' 


Fig. 


15 


Capak et al. (2004) 


z 


~ 5 


616 


z' 




25.8 


68 


R-i', i'-z' 


Fig. 


16, 20 


Ouchi et al. (2004) 


z 


~ 6 


12t 


z' 




29.5 


506 


i — z, V — z 


Fig. 


21 


Bouwens et al. (2005) 




~ () 


7(i7 






25.4 


12 




Fig. 


21 


Sliiuijisaku ot al. (2()()r)) 



'I' Bouwens et al. (2005) incorporated the faint data from the Hubble Ultra Deep Field (12 arcmin^) with two other data sets of 
shallower depths: the two GOODS fields (~ 170 arcmin^ each) and the two Ultra Deep Field-Parallel fields (20 arcmin^ in total). 
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Fig. 14. — Observed number counts of LBGs at 
2; ~ 4 uncorrected for incompleteness or contami- 
nation. The filled circles represent the BRi'-LBGs 
in this study. Comparison is made with the previ- 
ous measurements obtained by Ouchi et al. (2004) 
(open circles), Steidel et al. (1999) (open squares), 
Sawicki & Thompson (2006) (open pentagons), 
and Capak et al. (2004) (crosses). The error bars 
for this study, Ouchi et al. (2004), and Capak et 
al. (2004) reflect Poisson errors. The error bars 
for Steidel et al. (1999) and Sawicki & Thompson 
(2006) include an estimate of cosmic variance from 
their multiple fields as well as Poisson errors. 

aperture correction of —0.2 mag. The correction 
value, however, actually can range to as much 
as —0.4 mag for bright LBGs. A shift of the 
magnitudes of Ouchi et al. (2004) by —0.2 mag 
(= —0.4— (—0.2)) can settle the disagreement be- 
tween Ouchi et al.'s (2004) and our counts. The 
counts of Capak et al. (2004) are lower than our 
counts for 2 ~ 4 widely at faint magnitudes and 
for z ~ 5 at the faintest magnitude of Capak et al. 
(2004). The lower detection completeness of their 
data might be a major cause of this inconsistency; 
the 5(7 limiting magnitude of their data is about 1 
mag brighter than that of the SDF Project data. 
For 2; ~ 4 LBGs, our measurements agree with 
those of Steidel et al. (1999), although Steidel et 
al.'s (1999) measurements are limited to objects 
brighter than 25 mag. At faint magnitudes, how- 
ever, our measurements are higher than those of 
Sawicki & Thompson (2006). For 2 ~ 5 LBGs, 
over a wide range of bright magnitudes (m < 25), 
a discrepancy is found between the counts of Iwata 



Fig. 15. - Observed number counts of V-dropout 
LBGs at z ~ 5 uncorrected for incompleteness 
or contamination. The filled circles represent 
the Vi' z'-LiBGs in this study. Comparison is 
made with the previous measurements obtained 
by Ouchi et al. (2004) (open circles), Iwata et al. 
(2003) (open triangles), and Capak et al. (2004) 
(crosses). The error bars reflect Poisson errors. 

et al. (2003) and either ours or those of Capak et 
al. (2004) who surveyed almost the same field as 
Iwata et al. (2003); Iwata et al.'s (2003) counts 
are higher about a factor of 3 - 10. The selec- 
tion criteria of Iwata et al. (2003) are somewhat 
looser than ours. Therefore, as Ouchi et al. (2004) 
and Capak et al. (2004) claim, the LBG sample of 
Iwata et al. (2003) could be strongly contaminated 
by foreground galaxies. 

3.3. Completeness 

The completeness at a given apparent magni- 
tude and redshift is defined as the ratio of LBGs 
which are detected and also pass the selection cri- 
teria, to all the LBGs with the given apparent 
magnitude and redshift actually present in the 
universe. In general, completeness decreases as 
apparent magnitude goes fainter and as redshift 
goes away from the central redshift of the sam- 
ple. The completeness of each LBG sample is es- 
timated as functions of apparent magnitude and 
redshift through a Monte Carlo simulation. 

In the simulation, we generate artificial LBGs 
over an apparent magnitude range of 23.1 < m^/ < 
26.6 for BRi'-LBGs and 23.8 < m^/ < 25.8 for 
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Fig. 16. — Observed number counts of i?-dropout 
LBGs at 2; ~ 5 uncorrected for incompleteness 
or contamination. The filled circles represent 
the ffi'z'-LBGs in this study. Comparison is 
made with the previous measurements obtained 
by Ouchi et al. (2004) (open circles). The error 
bars reflect Poisson errors. 

yi'^'-LBGs and Ri'z'-LBGs with an interval of 
Am = 0.5, and over a redshift range of 2; = 3.0 < 
z < 5.9 with an interval of Az = 0.1. All LBGs arc 
assumed to have the same intrinsic spectrum de- 
scribed in §3.1, but with values oiE{B-V) = 0.0, 
0.1, 0.2, 0.3, and 0.4. Thus, 1, 200 (8 x 30 x 5) ob- 
jects for BRi'-hBGs and 750 (5 x 30 x 5) objects 
for Vi'z'-LBGs and Ri' z'-LBGs are generated. We 
assume that the shape of LBGs is PSF-likc and 
assign apparent sizes to the artificial LBGs so 
that their size distribution measured by SExtrac- 
tor (the peak FWHM is w l."00 - l."04) matches 
that of the observed LBG candidates. The arti- 
ficial LBGs are then distributed randomly on the 
original images after adding Poisson noise accord- 
ing to their magnitudes, and object detection and 
photometry are performed in the same manner as 
done for real objects. A sequence of these pro- 
cesses is repeated 100 times to obtain statistically 
accurate values of completeness. In the simula- 
tion, the completeness for a given apparent mag- 
nitude, redshift, and E{B — V) value can be de- 
fined as the ratio in number of the simulated LBGs 
which are detected and also satisfy the selection 
criteria, to all the simulated objects with the given 
magnitude, redshift, and E{B — V) value. We cal- 
culate the completeness of the LBG samples by 



taking a weighted average of the completeness for 
each of the five E{B — V) values. The weight is 
taken using the E{B — V) distribution function of 
2; ~ 4 LBGs derived by Ouchi et al. (2004) (open 
histogram in the bottom panel of their Fig. 20), 
which is corrected for incompleteness due to selec- 
tion biases. We plot the recovered artificial objects 
in the two-color diagrams by taking the number of 
objects with each E{B — V) value according to the 
adopted E{B — V) distribution fimction. We ver- 
ify that the distribution well resemble that of real 
objects. Hence, our modeling is reasonably realis- 
tic. The resulting completeness, p{m,z), is shown 
in Figure 17. 

In order to examine the uniformity of complete- 
ness over our images, we split the images into two 
regions (a center region and a periphery region), 
each with the same area, and recalculate the com- 
pleteness for each region. The non-uniformity is 
found to be small; the variation between the cen- 
ter region and the periphery region is within 10 
%. 

For each of the three LBG samples, magnitude- 
weighted redshift distribution functions are de- 
rived from p{m, z) (thick black lines in Figure 17) 
by averaging the magnitude-dependent complete- 
ness weighted by the number of LBGs in each 
magnitude bin. The average redshift, z, and 
its standard deviation, S2, are calculated to be 
z = 4.0 and = 0.3 for BRi'-hBGs, z = 4.7 and 

= 0.3 for T/i'z'-LBGs, and z = 4.9 and = 0.2 
for i?'i'2;'-LBGs. We also compute magnitude- 
weighted redshift distribution functions where the 
weighting is done by the number of spectroscopi- 
cally identified LBGs (dotted lines in Figure 11), 
to obtain the average redshift of Zgpoc = 4.0 
for BRi'-LBGs, z^pcc = 4.8 for yi'z'-LBGs, and 
Zspcc = 4.9 for Ri'z'-LBGs. These distribution 
functions are consistent with the redshift distribu- 
tions of spectroscopic samples shown in Figure 11. 

Although the above modeling well reproduces 
the observed distribution of objects in the two- 
color diagrams and the observed redshift distribu- 
tion functions, we further explore to what extent 
the results would be affected by adopting differ- 
ent models. To begin with, we recalculate the 
completeness using model spectra with two other 
ages, 0.01 Gyr and 0.5 Gyr, to find that the com- 
pleteness changes only a few percent. Then, we 
examine two extreme cases of dust extinction that 
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all model spectra have E{B — V) — 0.0 and 0.4. 
Changes in completeness are found to be negli- 
gibly small in either case except for Ri'z'-LBGs. 
Next, wc examine the effect of changing the ab- 
sorption by the intergalactic medium. If we shift 
the amount of attenuation to ilcr from the av- 
erage amount given in Madau (1995), we obtain 
quite different completeness values. However, the 
redshift distribution functions derived are unreal- 
istic, because they are strongly inconsistent with 
the distributions of spectroscopic objects. We also 
adopt Meiksin (2006) 's prescription for absorption 
by the intergalactic medium, to find completeness 
values similar to those based on Madau (1995)'s. 

3.4. Contamination by Interlopers 

We estimate the fraction of low-redshift inter- 
lopers in the LEG samples by a Monte Carlo sim- 
ulation as follows. For the boundary redshift, zq, 
between interlopers and LBGs, Zq = 3.5 is adopted 
for BRi'-LBGs, zq = 4.0 for Vi'z'-LBGs, and 
zo = 4.5 for Ri'z'-LBGs. 

We use objects in the Hubble Deep Field North 
(HDFN), for which best-fit spectra and photomet- 
ric redshifts are given by Furusawa et al. (2000), 
as a template of the color, magnitude, and redshift 
distribution of foreground galaxies, and generate 
929 artificial objects which mimic the HDFN ob- 
jects. The apparent sizes of the artificial objects 
are adjusted so that the size distribution recov- 
ered from the simulation is similar to that of the 
real objects in our catalogs. We distribute the 
artificial objects randomly on the original images 
after adding Poisson noise according to their mag- 
nitudes, and perform object detection and pho- 
tometry in the same manner as employed for real 
objects. A sequence of these processes is repeated 
100 times. In the simulation, the number of inter- 
lopers can be defined as the number of the sim- 
ulated objects with low redshift {z < zq) which 
are detected and also satisfy the selection criteria 
for LBGs. The number of interlopers expected in 
an LEG sample can then be calculated by multi- 
plying the raw number by a scaling factor which 
corresponds to the ratio of the area of the SDF 
(875 arcmin^) to the area of the HDFN multiplied 
by the repeated times (100 x 3.92 arcmin^). Fig- 
ure 18 shows the fraction of interlopers for our 
LEG samples as a function of magnitude. For the 
BRi'-LBG and Ri'z'-LBG samples, the fraction 



is found to be less than 5% at any magnitude. 
For the ^i'z'-LEG sample, the contamination is 
higher but at most ~ 20%. Most of the interlop- 
ers arc in the redshift range of 0.2 < z < 0.8, as 
predicted from Figures 4-6. The rest are objects 
at redshifts which are close to the boundary red- 
shift of each sample. The fraction of interlopers 
only at low redshifts (i.e., not at near the bound- 
ary redshifts) for each whole sample is around 1% 
for the BRi'-LBG and ffi'z'-LBG samples, and 
around 9% for the Vi'z'-LBG sample. 

The number density and the redshift distribu- 
tion of galaxies in the HDFN may be largely diS'er- 
ent from the cosmic averages because the HDFN 
is a very small field. However, the contamination 
by interlopers for each of the three LEG samples 
calculated above is very low. We therefore ex- 
pect that the uncertainty in contamination due 
to a possible (large) cosmic variance in the HDFN 
galaxies will not be a significant source of the er- 
ror in the luminosity functions of LEGs derived in 
the next section. 
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Fig. 17. — Completeness against the redshift of 
objects with different apparent magnitude for our 
LBG samples. In the panel of BRi'-LBGs, the 
red, magenta, green, cyan, blue, red, magenta, and 
green (from top to bottom) hues denote the com- 
pleteness for i' = 23.45, 23.95, 24.45, 24.95, 25.45, 
25.95, 26.45, and 26.95, respectively. In the panels 
of Vi' z'-LBGs and Ri' z'-LBGs, the red, magenta, 
green, cyan, and blue lines denote the complete- 
ness for z' = 24.15, 24.65, 25.15, 25.65, and 26.15, 
respectively. The thick black lines in all panels 
indicate the magnitude-weighted completeness. 
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Fig. 18. — Fraction of interlopers as a function 
magnitude for our LBG samples. 
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4. LUMINOSITY FUNCTIONS AT REST- 
FRAME ULTRAVIOLET WAVELENGTHS 

Wc derive the luminosity functions (LFs) of 
LBGs at z ~ 4 5 by applying the "effective vol- 
ume" method (Steidel et al. 1999). The Monte 
Carlo simulation in §3.3 yields the effective survey 
volume as a function of apparent magnitude, 



Vesim) 



f 

•I za 



P{'m,z)^^^dz, (4) 



where p{m, z) is the probability that a galaxy of 
apparent magnitude m at redshift z is detected 
and passes the selection criteria (i.e., the com- 
pleteness in §3.3 for magnitude m and redshift z), 

^^C^ is the differential comoving volume at red- 
ds: ^ 

shift z for a solid angle of the surveyed area (875 
arcmin^), and zq is the boundary redshift between 
low-redshift galaxies and LBGs. The effective vol- 
ume obtained for each magnitude bin is provided 
in Tables 3-5. 

The number density of LBGs corrected for in- 
completeness and contamination is computed as: 

, . s -/Vraw(TO) - iVintorloper(m) 

^^"^^ = ' 

where Nrawim) and iVinterioper(TO) are the number 
of LBGs detected and the number of interlopers 
which is estimated from the simulations, respec- 
tively, in an apparent magnitude bin of m. Tables 

3-5 also give iVraw(m) and A^interloper(m). 

The number densities of LBGs in faint bins, es- 
pecially the faintest bins, could be largely over- 
estimated owing to Eddington bias from objects 
in fainter bins and/or beyond the limiting magni- 
tude of the sample which have very large photo- 
metric errors and probably have a larger number 
density. We evaluate the effect of this bias for 
each LBG sample as follows. First, we construct 
a deeper sample of LBGs down to a magnitude 
limit fainter by 0.5 mag than the original value 
(for example, the new limit is i' < 27.35 for BRi'- 
LBGs). Then, using this deeper sample, we carry 
out Monte Carlo simulations similar to those made 
in §3.3, and derive Vesijn) to calculate (j)(rn) to a 
deeper magnitude. Prom the simulations, we com- 
pute for each magnitude bin the fractions of LBGs 
which scatter out to other magnitude bins. Fi- 
nally, we recalculate (j){m) by correcting the origi- 
nal number densities for the scatter along the mag- 



nitude axis. The number density of LBGs derived 
in this way is found to be consistent with the orig- 
inal calculation within statistical errors, down to 
the original limiting magnitude. This means that 
Eddington bias is negligibly small for our samples. 

The LF in the rest-frame ultraviolet (~ 1500 A) 
absolute magnitude, A'/uv, if^ obtained by convert- 
ing (j){in) into 4'{M-[jY{m)). For each LBG sample, 
we calculate the absolute magnitude of LBGs from 
their apparent magnitude using the average red- 
shift of the sample by 

Muv = m-h2.51og(l z) - 51og(dz,(2)) + 5 

-l-(TOtJV-TO), (6) 

where (II is the luminosity distance in units of 
pc. For apparent magnitude, m, we use i' magni- 
tude for _Bi?2'-LBGs, and z' magnitude for Vi' z'- 
LBGs and i?z'z'-LBGs. The last term of Eq.6, 
(wuv — "i)) is the difference between the magni- 
tude at rest-frame 1500 A and the magnitude in the 
bandpass being considered in the rest-frame. Us- 
ing the representative galaxy models given in §3.1, 
we find (muv — m) to be negligible, or no larger 
than 0.03 mag, over the redshift ranges selected by 
the color selection criteria we adopt. Thus we set 
(muv — m) to zero in Eq. 6. We assume that all 
of the LBGs have the average redshift of the sam- 
ple. Indeed, for each of the three LBG samples, 
varying the redshift of an object over the standard 
deviation of the redshift distribution of the sample 
changes its absolute magnitude by not more than 
0.1 mag. 

In Figures 19 and 20, the LFs at z ~ 4 and 
z ~ 5 are plotted respectively, along with those 
of other authors (See Table 2.). Note that Ca- 
pak et al. (2004) do not derive the LFs. Our LFs 
are in excellent agreement with those of Ouchi 
et al. (2004) up to their faintest magnitudes for 
both redshifts. This agreement appears reason- 
able, since Ouchi et al. (2004) and this study are 
based on the same field (the SDF) and adopt al- 
most the same selection criteria. Our LFs reach 
fainter magnitudes than theirs. However, there 
is a discrepancy between our z ^ 4 LF and that 
of Sawicki & Thompson (2006) in the Keck Deep 
Fields at the faint end; the faint-end slope of our 
LF is steeper than that of Sawicki & Thompson 
(2006). Gabasch et al. (2004) also derived the 
galaxy LFs at rest-frame ultraviolet wavelengths 
in the FORS Deep Field using photometric red- 
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Table 3 

Number of LBGs detected, number of interlopers, and effective survey volume for 

BRi'-LBGs. 





Magnitude Range {i') 


AT a 
-''raw 


AT- , b 
* interloper 




22.85 


- 23.35 


3 





1.74 X 10^ 


23.35 


- 23.85 


12 





1.63 X 10<^ 


23.85 


- 24.35 


68 





1.49 X 10^ 


24.35 


- 24.85 


231 





1.39 X 10^ 


24.85 


- 25.35 


447 





1.22 X 10*5 


25.35 


- 25.85 


763 


8.9 


1.15 X 10*^ 


25.85 


26.35 


1093 


24.6 


8.77 X 10^ 


26.35 


- 26.85 


1191 


26.8 


5.70 X 10^ 



* Number of LBGs detected. 

^ Number of interlopers estimated from our simulations (see §3.4). 
Effective survey volume for 875 arcmin^ in units of Mpc^. 



Table 4 

Number of LBGs detected, number of interlopers, and effective survey volume for 

yi'^'-LBGs. 





Magnitude Range {z') 


N ^ 


^ interloper 




23.55 


~ 24.05 


1 





1.34 X 10*^ 


24.05 


- 24.55 


12 





1.31 X 10*^ 


24.55 


- 25.05 


56 





1.31 X 10'^ 


25.05 


- 25.55 


152 


29.0 


1.17 X 10^ 


25.55 


- 26.05 


318 


22.3 


9.66 X 10^ 



Number of LBGs detected. 
^ Number of interlopers estimated from our simulations (see §3.4). 
Effective survey volume for 875 arcmin^ in units of Mpc^. 
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Table 5 



Number of LBGs detected, number of interlopers, and effective survey volume for 

Ri'z'-LBGs. 





Magnitude Range {z') 




interloper 




24.05 


- 24.55 


3 





3.32 X 10'^ 


24.55 


- 25.05 


18 





4.60 X 10^ 


25.05 


- 25.55 


66 


2.2 


5.54 X 10^ 


25.55 


- 26.05 


153 


2.2 


5.69 X 10^ 



Number of LBGs detected. 

Number of interlopers estimated from our simulations (see §3.4). 
^ Effective survey volume for 875 arcmin^ in units of Mpc^. 



shifts and obtained a similar result to that of Saw- 
icki & Thompson (2006). The reason for this in- 
consistency is not clear. We note, however, that 
our raw counts before the correction of incomplete- 
ness already exceed the Sawicki & Thompson's 
(2006) corrected counts, and thus that our higher 
values are not due to overcorrections. Addition- 
ally, our data have a sky coverage five times larger 
than that of the total of the five Keck Deep Fields. 
Thus, our LF is more robust against a possible cos- 
mic variance on large scales, although the five sep- 
arate Keck Deep Fields are less affected by cosmic 
variance than a single continuous field of the same 
sky coverage. At 2 ~ 5, our LF is different from 
that of Iwata et al. (2003) at bright magnitudes. 
One possible explanation for this difference is that 
Iwata et al. (2003) may select a large number of 
contaminants as mentioned in §3.2, resulting in 
the higher number density of bright LBGs. 

We find a good agreement between our 2 ~ 5 
LFs obtained from the Vi'z'-LQG sample and the 
Ri' z'-LQG sample. The redshift ranges of these 
two samples are almost the same. This consistency 
in the results from the two independent man- 
ners provides strong support for LFs we obtained. 
As the Vi' z'-IjQG sample contains a statistically 
larger number of objects and has higher complete- 
ness, we only use results from the V^i'z'-LBG sam- 
ple in the following analysis for z ~ 5. These LFs 
obtained at 2; ~ 4 and 5 are found to be well re- 
produced by a semi-analytic model combined with 
high-resolution N-body simulations (Nagashima et 



al. 2005) when the observational selection effects 
are taken into account. We should, however, point 
out that the Lyman-break technique cannot se- 
lect all star-forming galaxies at targeted redshifts. 
On the basis of spectroscopy of galaxies in an 
optically-selected, flux-limited sample, Le Fevre et 
al. (2005) claim that a large fraction of galax- 
ies at 2; ~ 3 - 4 lie outside the selection bound- 
aries for LBGs in two-color diagrams. In addition, 
submm-selected galaxies (e.g., Small, Ivison, & 
Blain 1997) and NIR-selected star-forming galax- 
ies (e.g., Frarrx et al. 2003) are often too faint in 
optical wavelengths to be selected as LBGs, be- 
cause of strong absorption by dust. Such dusty 
star-forming galaxies could contribute to the cos- 
mic SFR density as much as do LBGs (e.g., Chap- 
man et al. 2005). Therefore, the far UV luminos- 
ity function and the star formation rate density 
derived in this paper are for galaxies with modest 
dust extinction and passing selection criteria for 
LBGs, and thus they may be greatly modified if 
all star-forming galaxies are included. 

We fit the LFs with the Schechter function 
(Schechter 1976): 

<j){L)dL 



L* 



or expressed in terms of absolute magnitude: 

(l){Muy)dMuv 
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Fig. 19. — UV luminosity functions (LFs) of LBGs 

at z ^ 4. Our data arc shown by the fiUcd circles. 
The open circles, open squares, and crosses are 
from Ouchi et al. (2004), Steidel et al. (1999), and 
Sawicki & Thompson (2006), respectively. The 
error bars for our data and Ouchi et al.'s (2004) 
reflect Poisson errors, and those for the other two 
data include both Poisson errors and an estimate 
of field-to-field variance from their multiple fields. 



X exp (-l0-0-4(^uv-M5v) j rfMuv, (8) 

where (f)*, L* (Af{jy), and a are parameters to be 
determined from the data. The parameter cj)* is 
a normalization factor which has a dimension of 
number density, L* is a "characteristic luminos- 
ity" (with an equivalent "characteristic absolute 
magnitude", M^y), and a gives the slope of the 
luminosity function at the faint end. For the 2; ~ 5 
LF, whose a has a rather large uncertainty, we 
evaluate </>* and M^y with a fixed to the best- 
fit value for the ^; ~ 4 LF, besides fitting with all 
three parameters left free. The best-fit parameters 
obtained are listed in Table 6. 

5. EVOLUTION OF THE COSMIC STAR 
FORMATION ACTIVITY OVER < 

z<6 

5.1. Evolution of the Luminosity F\inction 

We plot the LFs at 2; ~ 4 and 5 obtained in 
this study with those of LBGs at z 3 (Sawicki 
& Thompson 2006) and z ~ 6 (Bouwens et al. 



Fig. 20.— UV luminosity functions (LFs) of LBGs 

at z ^ 5. Our data from Vi' z'-LiBG sample and 
Ri'z'-LBG sample are shown by the filled circles 
and filled triangles, respectively. The open circles 
and open triangles are from Ouchi et al.'s (2004) 
Vi'z'-LBG sample and Ri'z'-LBG sample, respec- 
tively. The crosses are from Iwata et al. (2003). 
The error bars reflect Poisson errors. 

2004a; Shimasaku et al. 2005 from observation 
of the SDF in new filters sensitive to ~ 1 /xm) 
in Figure 21. The LFs of UV (1500A)-selected 
galaxies at z (Wyder et al. 2005) and z ~ 
1 (Arnouts et al. 2005) are presented as well for 
reference. Figure 21 reveals clear evolution of the 
LF beyond z ~ 4. From z ~ 4 to z ~ 3, we find 
no significant evolution over the whole luminosity 
range. Sawicki & Thompson (2006) and Gabasch 
et al. (2004) detect an increase in faint galaxies 
across this redshift interval, but such an increase is 
not seen in our data. At lower redshifts, luminous 
galaxies strongly decrease in number from z > 2 
to z ~ 0, as has been found by previous studies. 

It is unlikely that the evolution in LF seen in 
Figure 21 is caused by the differences in selection 
criteria among the different samples. We find, us- 
ing spectra of model galaxies, that the criteria for 
LBGs at z ~ 3, 4, and 5 select galaxies over simi- 
lar ranges of age and E{B — V). Since i'-dropout 
galaxies at z 6 arc selected with i' — z' color 
alone, they can have a wide range of E{B — V). 
However, since no observation has reported a sig- 
nificant increase in dust extinction at z > 5, it 
is reasonable to assume that the i'-dropout galax- 
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Table 6 

Luminosity function parameters 



Sfimplo NajiK" : .s^ o* "Mix: '] -^-^i v " 



BRi'-LBGs 


4.0 


0.3 




-3 


-21 14+014 


-1 00+0.09 
■'■•°'^-0.09 


Vi'z'-LBGs 
Vi'z'-LBGs 


4.7 
4.7 


0.3 
0.3 


0.58+J-^^ X 10 
1.23tg;^^ X 10- 


-3 
-3 


-21.09t0-?4 
-20.72l°i| 


9 O1+0.68 
^•^-•■-0.60 

-1.82 (fixed) 



ies are a similar population to LBGs at lower red- 
shifts. Galaxies at z < 1 are selected in completely 
different ways from the Lyman-break technique. 
However, the differences in LF between z > 3 and 
2: < 1 seen in Figure 21 appear to be too large 
to be accounted for by possible selection effects. 
The dimming of LF seen for z <1 will at least be 
qualitatively correct. 

Figure 22 shows the evolution of Schechter pa- 
rameters over < z < 6. The filled circles repre- 
sent our measurements. The other symbols show 
measurements taken from the literature. Note 
that the measurements at z 2 and z ^ 3 from 
Arnouts et al. (2005) may have large uncertainties, 
because they are based solely on the HDFs which 
cover only a tiny area. The normalization factor, 
(/>*, and the faint-end slope, a, appear to be almost 
constant over the redshift range of 0.5 < z < 6, 
followed by a significant change at lower redshifts. 
We note that a at z ~ 5 is fixed to the value at 
.2 ~ 4 in evaluating the other two parameters but 
the value at z ~ 4 is similar to that at z ~ 6. From 
2:~0.5toz~0, 0* rises by a half order of magni- 
tude and a becomes shallower from —1.6 to —1.2. 
On the other hand, the characteristic magnitude, 
M*, strongly and non-monotonically evolves over 
< z < 6. It brightens by about 1 mag from 
z ~ 6 to z 4, remains unchanged to z 3, and 
then a significant fading of about 3 mag occurs 
from 2 ~ 3 to z ~ 0. In other words, the evolution 
of the LF seen in Figure 21 is accounted for by 
a change in M* with cosmic time. We point out, 
however, that there are different results about the 
LFs Sit z ^ 4 and z 5, as mentioned in §4. Some 
groups claim that the faint-end slope is shallower 
at higher redshifts (Iwata et al. 2003; Gabasch et 
al. 2004; Sawicki & Thompson 2006). The charac- 
teristic magnitude of Sawicki & Thompson (2006) 



at z ^ 4 is consistent with ours. Iwata et al. (2003) 
obtained an about 1 mag brighter value at z ~ 5 
than ours. 

5.2. Cosmic Star Formation Rate Density 

5.2.1. Cosmic Star Formation Rate Density be- 
yond z ^ 3 

The ultraviolet luminosity density, pi,, at 2 ^ 4 
and 5 is derived by integrating the LFs obtained 

in §4 as: 

POO 

PL = L^{L)dL. (9) 

We consider two kinds of luminosity density. One 
is the "observed" luminosity density, p'^^, which 
is derived from integration down to the faintest lu- 
minosity of the sample, Lmin = ^min- This quan- 
tity, p'^^, gives a lower limit of the true luminos- 
ity density. The faintest absolute magnitude of 
the sample at each redshift is Muv — —19.2 for 
z ~ 4 (BRi'-LBGs) and Muv = -20.3 for z ~ 5 
(Fi'z'-LBGs). The other is the luminosity density 
derived by integrating the LFs down to the same 
luminosity limit for all the measures. This is im- 
portant to properly account for possible evolution. 
Following Steidel et al. (1999), we choose the limit 
to be imin — 0.1L*^3, where Ll^^ is the charac- 
teristic luminosity of the LF of z ~ 3 LBGs given 
by Steidel et al. (1999). The luminosity O.lLJ^g 
corresponds to Afuv = —18.6. We refer to this 
quantity as p^>°-^^^=3 hereafter. When calculat- 
ing p^^*^'^^"""^, we extrapolate the Schechter func- 
tions derived in §4 to fainter magnitudes than the 
observation limits. This extrapolation is modest 
for the z ~ 4 LBG sample (i?i?i'-LBGs), since 
L^Sn of this sample (Muv = -19.2) is close to 
0.1L*^3, On the other hand, the extrapolation is 
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Fig. 21. — UV luminosity functions (LFs) of LBGs 
at z 3 (blue open circles: Sawicki & Thompson 
2006), 

z ^ 4 (green filled circles: this study [BRi'- 
LBGs]), z ^ 5 (red filled circles: this study [Vi'z'- 
LBGs]), and z ~ 6 (cyan crosses: Bouwens et al. 
2005; cyan filled triangle: Shimasaku et al. 2005). 
The error bars reflect Poisson errors. The best-fit 
Schechter function is also shown with the solid line 
for each data set except for that of Shimasaku et 
al. (2005). The black dashed line and the black 
dotted line indicate the LFs of UV-selected galax- 
ies at z ~ (Wyder et al. 2005), and z ~ 1 
(Arnouts et al. 2005), respectively. 

large for the z ~ 5 LBG sample {Vi' z'-LBGs); 
MjjY — —18.6 is fainter than the limiting magni- 
tude of this sample by 1.7 mag. Table 7 presents 
pt^ and 

The ultraviolet luminosity density can be used 
to measure the star formation rate (SFR) density 
in the universe. We compute the cosmic SFR den- 
sity using the relationship between SFR and ultra- 
violet luminosity, Luv, given by Kennicutt (1998): 

SFR [Mq yr-i] 

= 1.4 X IQ-^^LuY [ergs s'^ Hz'^]. (10) 

This formula assumes a Salpeter initial mass func- 
tion with 0.1 Mq < M < 100 Mq. The SFR 
density is corrected for dust extinction follow- 
ing Hopkins (2004), who used the dust extinc- 
tion formula by Calzetti et al. (2000) and assumed 
E{B - y) = 0.128 at all redshifts. 

In Figure 23, we show the cosmic SFR densities 
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Fig. 22. — Evolution of the luminosity function 
parameters with time. The filled circles repre- 
sent the measurements from this study. The other 
symbols show measurements taken from the lit- 
erature: open squares from Wyder et al. (2005), 
filled squares from Arnouts et al. (2005; GALEX), 
filled triangles from Arnouts et al. (2005; HDF), 
open circles from Sawicki & Thompson (2006), and 
crosses from Bouwens et al. (2005). For z ~ 5, we 
fixed a to the value at z ~ 4 (q = —1.82) in eval- 
uating the other two parameters. 



obtained in this study as a function of redshift 
(filled circles), comparing with those calculated 
from LFs of LBGs at 3 < z < 6 in the literature 
(other kinds of symbols). We similarly integrate 
their LFs down to their observed limiting lumi- 
nosities to obtain pf^^ and down to L = 0.1L*^3 
to obtain p^'^'^'^^''=^ , convert them to the SFR 
densities using Eq.(lO), and correct for the same 
amount of dust extinction to provide a fair com- 
parison. The black symbols indicate Pj^'^'^^ 
while the cyan symbols are for the lower limit, 
p'j^^- Note that the limiting luminosity of the LF 



at z ~ 6 by Bouwens et al. (2005) is fainter than 
0.1L*^3, resulting in pf'^ > p^^° "'^'^'=\ It can 

be seen that p^^*'^ is almost constant over 
3 ^ 2 ^ 5, and perhaps turns to decrease at a 
certain point between 2^5 and 6. It should 
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Table 7 
uv luminosity densities 



PL [ergs~^ Hz~^ Mpc~^] 



Sample Name z Sz Pl 

(Muv < Miimit* ) (Muv < -18.6) 



Bi?z'-LBGs 4.0 0.3 2.2lt;^-^i x lO^^^ '^■^'^-I'll x 10^' 

r/'-'-LBGs 4.7 ().:-! 4.:iSi;j;;]'^ X 1.85+;;;';;;^ x lO^' 



* Observed limiting magnitudes; Muv = -19.2 for BRi' -IBGs and Muv = -20.3 
for Vi'z'-hBGs. 



be emphasized that the measurements of the SFR 
density at 2; ~ 4 and 2; ~ 5 are largely improved 
by this study. At z ~ 4, the observed SFR den- 
sity is very close to p^^'^^'"^, since our BRi'- 
LBG sample reaches very close to 0.1L*^3. Thus, 
there is little uncertainty in our measurement of 
^L>o.iL^^3 g^i; 2 '-vj 4, while the previous measure- 
ments of SFR density integrated to 0.1L*^3 have 
uncertainties of a factor of two or more due to a 
large extrapolation of the LF. From this figure, we 
can robustly conclude that the cosmic SFR den- 
sity does not drop from z ~ 3 to z ~ 4, since 
the "observed" SFR density at z 4 is not lower 
than Pj^^'^ '"^ at z ~ 3. At z ~ 5, large errors 
owing to the uncertainty in the shape of the LFs 
at faint magnitudes still remain. However, from 
our observed SFR density, we can put a stronger 
constraint on the "decrease" in the cosmic SFR 
density from z ~ 4 to z ^ 5; the decrease, if any, 
cannot be larger than a factor of about five. 

5.2.2. Luminosity Dependent Evolution of the 
Cosmic Star Formation Rate Density 
over < z < 6 

We examine the evolution of the cosmic SFR 
densities contributed from galaxies with different 
magnitudes over the redshift range of < z < 6. 
Figure 24 compares their evolutionary behavior. 
The black, red. magenta, green, and blue symbols 
indicate the cosmic SFR densities from galaxies 
with magnitudes in the range M < M*^3 — 1.5, 

m;^3-i.5 < m < m;^3-o.5, m;^3-o.5 < m < 

M;^3 + 0.5, M;^3 + 0.5 < M < M;^3 + 1.5, and 



M*^3 -I- 1.5 < M < M*^3 -t- 2.5, respectively. The 
filled circles represent our measurements, and the 
other symbols show the values similarly calculated 
using the LFs taken from the literature. There is a 
remarkable dependence of the evolutionary behav- 
ior of the SFR density on luminosity, as is previ- 
ously noticed by Shimasaku et al. (2005) who com- 
pare the evolution of the total cosmic SFR den- 
sity and the SFR density from galaxies brighter 
than M = -21.3 to find that the SFR density 
from bright galaxies drastically changes with time. 
From Figure 24, it is shown that the SFR den- 
sity from fainter galaxies evolves more mildly. In 
addition, the peak position varies with luminos- 
ity: the contribution from fainter galaxies peaks 
at earlier epochs. For example, the SFR density 
for galaxies with M*^3 - 1.5 < M < M*^3 - 0.5 
rises by about an order of magnitude from z ^ 6 
to z ~ 4, remains unchanged to z ~ 3, and then 
drops by more than three orders of magnitude to 
the present epoch, thus having a sharp peak at 
z ~ 3 - 4. On the other hand, the SFR density 
for galaxies with M*^3 -|- 1.5 < M < M;^3 -|- 2.5 
is almost constant within a factor of about two 
at z > 2 with a blunt peak at z ~ 4 - 5. Iwata 
et al. (2003), Gabasch et al. (2004), and Sawicki 
& Thompson (2006) have reported different evo- 
lutions of the LBG LF from z ~ 4 or 5 to 3 from 
what we find in the SDF. At 3 < z < 5, the lu- 
minosity density of these authors shows an oppo- 
site trend that fainter galaxies evolve more rapidly. 
However, it is at z < 3 and z > 5 that most drastic 
evolution occurs in the luminosity density. Thus, 
using the LFs from these authors does not quali- 
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Fig. 23. — Cosmic star formation rate (SFR) den- 
sity as a function of redshift. The black sym- 
bols indicate the SFR densities calculated by inte- 
grating the luminosity functions down to Muv = 
— 18.6 (equivalent to L = 0.1L*^3 from Stei- 
del et al. 1999), while the cyan symbols show 
the lower limits, i.e., the contribution from ac- 
tually observed galaxies only. The filled circles 
represent the measurements from this study; the 
other points come from Steidel et al. (1999) (open 
squares), Iwata et al. (2003) (open triangles), 
Ouchi et al. (2004) (open circles), Gabasch et 
al. (2004) (open hexagons), Sawicki & Thomp- 
son (2006) (open pentagons), and Bouwens et al. 
(2005) (crosses). 

tatively change the results on the overall evolution 
over < 2 < 6 obtained here. 

This luminosity-dependent evolution of the cos- 
mic SFR density reflects the evolution of the LF 
discussed in §5.1. The strong increase in the SFR 
density for galaxies brighter than M*^3 — 0.5 from 
z ~ 6 to z ~ 3 is due to the brightening of M* with 
time. The SFR densities for these bright galaxies 
then drop to the present epoch, as M* fades dur- 
ing the same period. The evolution of the SFR 



densities for galaxies fainter than Af*^3 — 0.5 is 
mild, since they are less sensitive to the change in 
M* . The earlier peak position of the SFR density 
of fainter galaxies suggests that faint galaxies are 
more dominant in terms of the cosmic SFR density 
at earlier epochs. 

We cannot completely rule out the possibility 
that the luminosity-dependent evolution in the 
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Fig. 24. — Cosmic star formation rate (SFR) den- 
sity as a function of redshift over the redshift 
range of < 2 < 6, demonstrating its luminos- 
ity dependent evolution. The black, red, magenta, 
green, and blue symbols indicate the SFR density 
from galaxies with M < M*^^, M*^^ - 1.5 < 
M < Af - 0.5, Af;^3 - 0.5 <M < M*^3 -I- 0.5, 
M* + 0.5 < M < M* + 1.5, and M* + 1.5 < 



z=3 ■ 

M < A'/*^3 -f 2.5, respectively. The filled circles 
represent the measurements from this study. The 
open squares comes from Wyder et al. (2005), the 
filled squares from Arnouts et al. (2005; GALEX), 
the filled triangles from Arnouts et al. (2005; 
HDF), the open circles from Sawicki & Thomp- 
son (2006), 

and the crosses from Bouwens et al. (2005). 



cosmic SFR density seen above is not real but 
reflects luminosity-dependent evolution of some 
other properties such as dust extinction. If, for 
instance, galaxies have the smallest dust extinc- 
tion at z 3 ~ 4, then the observed M* will be 
the brightest at these redshifts, as seen in Fig- 
ure 22, even if the intrinsic, dust-corrected M* 
does not change with redshift. In this case, the 
behavior of the luminosity-dependent SFR den- 
sity will be qualitatively similar to that found in 
Figure 24. However, no observation has reported 
strong evolution of E{B — V) at least for LBGs. 
The E{B — V) distribution of LBGs appears to 
be unchanged over z ~ 5 and 3 (Iwata et al. 2003; 
Ouchi et al. 2004). In addition, E{B-V) is known 
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to be independent of apparent M (Adelberger & 
Steidel 2000; Ouchi et al. 2004). 

5.3. Specific Star Formation Rate 
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We explore the evolution of what we define as 
the specific SFR here, i.e., the cosmic SFR per 
unit baryon mass in dark haloes. The specific 
SFR means the efficiency of star formation aver- 
aged over the all dark haloes present at a given 
redshift. The specific SFR is computed by divid- 
ing the cosmic SFR density, p^^^^'^^ , by the 
average mass density of baryons confined in dark 
haloes, pb- Here, pb is calculated by 

Pb{z) = / n{M,z)MdM, (11) 

where fi^ = 0.05 is the density parameter of 
baryons, n(M, z) is the mass function of dark 
haloes at redshift z predicted by the standard Cold 
Dark Matter model with erg = 0.9, M^in{z) is 
the halo mass corresponding to a virial temper- 
ature of 1 X 10"^ K. We assume here that gas does 
not cool (and thus stars are not formed) in dark 
haloes with T < 1 x 10"^ K. Note that Mmin(z) is 
much smaller than the estimated mass of haloes 
which host the faintest galaxies we account here, 
i.e., galaxies with OAL^^^. 

Figure 25 shows the evolution of the specific 
SFR over < z < 6. The filled circles represent 
our measurements and the other symbols show the 
values from the literature. The dotted line in the 
figure shows (1 + z)"^ evolution. We find that the 
specific SFR increases in proportion to (1 -|- z)^ 
up to z ^ 4. In other words, the star formation 
efficiency in the universe rises with redshift. At 
> 4, on the other hand, the specific SFR appears 
to decrease. 

Here, according to the standard Cold Dark 
Matter model, the average internal density in 
physical units of baryons in dark haloes virialized 
at a given redshift, (Pb^), is known to be approx- 
imately proportional to the average mass density 
of the universe at that time, and the average mass 
density of the universe changes as (l + z)^. There- 
fore, it implies that the specific SFR increases ap- 
proximately in proportion to (Pb^) up to z ^ 4. 
To put it another way, since the specific SFR is 
the internal density of SFR in dark haloes, (psfr), 
divided by (Pb^), this result can be expressed in 



Fig. 25. Specific star formation rate (SFR) ver- 
sus (1 + z) over < z < 6. Symbols are the same 
as in Figure 24. The dotted line shows evolution 
in power-law, specific SFR oc (1 -|- z)^. 

terms of (psp^) and (Pb^) as: 

{p^¥k) oc (p°H)^ (12) 

It is surprising that the star formation in galaxies 

obeys such a simple power law over 90% of the 
cosmic history {0 < z < 4). Equation (12) resem- 
bles the Schmidt law (Schmidt 1959) of nearby 
disk galaxies, a relationship between the disk- 
averaged SFR and cold gas surface densities of 
5^SFR = const. X with A/' ~ 1.4 (Kennicutt 
1998). 

We present here one possible explanation of the 
observed behavior of the specific SFR. The ana- 
lytical modeling of galaxy formation based on the 
CDM model by Hernquist & Springel (2003) shows 
that the gas cooling rate, Mcooi/-/Wdh, of dark 
haloes roughly scales with H{z)^, or equivalently 
with (1 -I- z)^ as H{z) is approximately propor- 
tional to (1 -I- z)^^^ in our cosmology. They also 
find that the star formation rate in dark haloes 
scales with the cooling rate at low and interme- 
diate redshifts, while it becomes to be no longer 
dependent on the cooling rate at redshifts larger 
than a certain limit, where the internal density of 
haloes is so high that the cooling time is shorter 
than the consumption rate of cooled gas into stars. 
If we adopt this modeling, then our result implies 
that the cosmic SFR is primarily determined by 
the gas cooling rate at z < 4, while it is governed 



25 



by the conversion rate of cooled gas into stars at 
larger redshifts. 

Of course, the cooling rate may not be the ma- 
jor factor determining the c;osmic star formation 
rate. There are many other factors on which the 
star formation in galaxies can largely depends, 
such as the feedback by galactic winds, the UV 
background radiation field, and environmental ef- 
fects like galaxy merging and interactions. Since 
the relative importance of these factors will vary 
with the mass of dark haloes, it will be interesting 
to measure the specific SFR as a function of halo 
mass. 

6. SUMMARY AND CONCLUSIONS 

We investigated the luminosity functions and 
star formation rates of Lyman-break galaxies 
(LBGs) at z ^ 4 and 5 based on the optical imag- 
ing data obtained in the Subaru Deep Field (SDF) 
Project. The SDF Project is a program conducted 
by Subaru Observatory to carry out a deep and 
wide galaxy survey in the SDF. 

Three samples of LBGs in a contiguous 875 
arcmin^ area are constructed. One consists of 
LBGs at 2 ~ 4 selected with the B — R vs R—i' 
diagram {BRi'-LQGs). The other two consist of 
LBGs at z 5 selected with two kinds of two- 
color diagrams: V — i' vs i' — z' (^i'z'-LBGs) and 
R-i' vs i'-z' {Ri'z' -LBGs). The number of LBGs 
detected is 3, 808 for Bi?i'-LBGs, 539 for Vi'z'- 
LBGs, and 240 for Ri'z'-hBGs. The adopted se- 
lection criteria are proved to be fairly reliable by 
the spectroscopic data. We use Monte Carlo simu- 
lations to estimate the completeness and the frac- 
tion of contamination by interlopers. The redshift 
distribution functions obtained by the simulations 
are well consistent with the redshift distributions 
of spectroscopically identified objects in the LBG 
samples. 

We derived the luminosity functions at rest- 
frame ultraviolet wavelengths down to as faint as 
Muv = -19.2 at z ~ 4 and Muv = -20.3 at 
z ~ 5. The cosmic star formation rate (SFR) den- 
sity was measured from the ultraviolet luminosity 
density derived by integrating the luminosity func- 
tions. On the basis of the observed SFR density 
and the standard Cold Dark Matter model, the 
cosmic SFR per unit baryon mass in dark haloes, 
i.e., the specific SFR, is computed. Combining 



the results obtained in this work with those taken 
from the literature, we reached the following con- 
clusions: 

(i) Clear evolution of the luminosity function is 
detected. The characteristic magnitude M* 
changes strongly over the redshift range of 
< z < 6. It brightens by about 1 mag from 
z ~ 6 to z ~ 4, remains unchanged to z ~ 3 
and then fades by about 3 mag from z ^ 3 to 
z ~ 0. On the other hand, the normalization 
factor, (j)* , and the faint-end slope, a, are 
almost constant over 0.5 < z < 6. 

(ii) The measurements of the SFR density at 
z 4 and z ~ 5 are greatly improved, since 
the luminosity functions are derived down 
to very faint magnitudes. Hence, the evolu- 
tionary behavior of the SFR density is more 
constrained. The cosmic SFR density does 
not drop from z ~ 3 to z ~ 4, and the de- 
crease from z 4 to z ~ 5, if any, cannot be 
larger than a factor of about five. 

(iii) The SFR density contributed from brighter 
galaxies changes more drastically with cos- 
mic time. The contribution from galaxies 
brighter than M*^^ — 0.5 has a sharp peak 
around z = 3 4, while that from galax- 
ies fainter than M*^3 — 0.5 evolves relatively 
mildly with a broad peak at earlier epoch. 
This Imuinosity-dependent evolution of the 
cosmic SFR density suggests that the galaxy 
population that contributes to the total SFR 
density in the universe varies with time, re- 
flecting the evolution of the luminosity func- 
tion. 

(iv) The specific SFR is found to scale with red- 
shift as (1 -|- z)^ up to z 4, implying that 
the efficiency of star formation is on aver- 
age higher at higher redshift in proportion to 
the cooling rate within dark haloes. It seems 
that this is not simply the case at z>A. 

We would like to thank the referee for valuable 
comments and suggestions. We are deeply grate- 
ful to the Subaru Telescope staff for their devoted 
support in this project. 
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